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ABSTRACT 



A chondrule formation theory is presented where the chondrule formation zone is located 
within 0.1 AU of the protosun. This hot, optically thick, inner zone of the solar accretion disk 
is coincident with the formation region of the protosolar jet. It is suggested, that chondrules 
are ablation droplets produced by the interaction of the jet wind with macroscopic bodies 
that stray into the jet formation region. 

Provided these droplets are small enough, they will be swept up by the jet wind and 
subsequently cool at approximately the same rate as the expanding gas in the jet. There is 
a critical gas density (~ 10 -11 g cm -3 ) below which chondrules will undergo large, damped 
oscillations in altitude and thereby suffer reheating. 

We claim that it is in the cooler, high altitude regions above the mid plane of the inner 
accretion disk that compound chondrules are formed, and the collisional fragmentation of 
chondrules take place. Since these processes take place in the jet flow, one can make a 
prediction for the expected structure of triple compound-chondrules. For such chondrules, 
it is suggested that two "relatively large" secondary chondrules will avoid each other. This 
prediction is valid only if the gas-flow is sufficiently laminar or if the "spin-down time" for 
a double compound chondrule is less than the inter-chondrule collision time. 

The model assumes that particles, ranging in diameter from 1 /zm to 1 cm, can be 
ejected from the inner-accretion disk by the jet flow, and that the angular momentum of 
this material is sufficient to eject it from the jet flow. Given these assumptions, any material 
so ejected, will fly across the face of the accretion disk at speeds greater than the escape 
velocity of the system. This material can only be recaptured through the action of gas 
drag. Such a capture process naturally produces aerodynamic size sorting of chondrules and 
chondrule fragments, while the ejection of refractory dust provides a possible explanation 
for the observed complementarity between matrix and chondrules. 

This transfer of material will result in the loss of angular momentum from the upper 
atmosphere of the outer accretion disk and thereby facilitate the accretion of matter onto 
the protosun. 

PREFACE 

In this paper we discuss a relatively new chondrule formation model, which uses a hypothetical 
protosolar jet as the main chondrule formation mechanism. As this model is in the earliest stages of 
development, this particular monograph should not be treated as the forever, enduring, last word on the 
subject. Instead, we hope that our discussion will demonstrate the explanative and predictive power of 
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this model, and possibly prompt our colleagues to view chondrule formation from a different point of 
view as compared to other theories. 

Some readers would argue that fewer, not more, chondrule formation theories are required. Cer- 
tainly with nearly twenty different published theories, how can one determine which theory, if any, is 
the correct one? The answer, of course, is to embrace the theory that can not only explain many of the 
observed data, but also predict new facts. So in this study, we will attempt to explain: 

(1) the chondrule formation process; 

(2) the chondrule cooling rate; 

(3) the collisional fragmentation of chondrules and the formation of adhering compound chondrules, 
where the secondary (and smaller) chondrule was plastic at the time of collision; 

(4) chondrule reheating; 

(5) the chondrule upper size limit; 

(6) chondrule size sorting, 
and 

(7) the complementary chemical composition of matrix and chondrules. 

We will also discuss a possible link between protostellar jets and disk accretion, where a protostellar 
jet may inject material with low angular momentum into the upper atmosphere of an accretion disk, 
thereby enhancing stratified accretional flow, i.e., the upper layers of the disk accrete onto the protostar 
more readily than do the layers adjacent to the mid plane of the accretion disk. 

As for predictions, we will describe the expected structure of triple compound chondrules, i.e., where 
two small secondaries are adhering to a large primary. We suggest that, for such a compound chondrule, 
the two secondary chondrules will tend to "avoid" each other. We give the mathematical expression for 
the minimum avoidance angle between the two secondaries and discuss under what circumstances this 
prediction will break down. We also show that the experimental data collected so far (Unfortunately, 
for only four such chondrules) have structures that are consistent with the predicted structure. It is our 
fond hope, that we will inspire some of our colleagues to collect much more data, so that the validity 
or otherwise of this prediction can be confidently determined. 



1. JET ABLATION 

In the late 1970's, it was found that protostellar systems formed bipolar outflows of material, i.e., 
the gas flows travelling in two opposing directions, approximately perpendicular to the plane of the 
accretion disk (for a review, see Beckwith and Sargent 1993). These flows were detected via the CO 
rotation lines and consisted mostly of molecular material. The flow speeds of the bipolar outflows were 
found to be around 20 km s -1 , and by dividing their length by the observed flow speed, one could 
deduce a "dynamic" lifetime for these flows of around 10 5 years. Later, unbiased surveys of protostellar 
systems suggests that all protostellar systems undergo some form of bipolar outflow stage (Fukui et al. 
1993). 



- 3- 



Often, but not always, one can find within a bipolar outflow a faster, more collimated flow known as 
a protostellar jet. Protostellar jets are usually detected in the SII, 01 and Ha lines and are consequently 
referred to as optical jets. They have observed wind speeds in the range of 100 - 400 km s" 1 , they eject 
a large amount of gas (total mass 10~ 3 — 10 _1 M ), are quite energetic (total kinetic energy 10 44 — 10 46 
erg), long-lived phenomena (10 6 — 10 7 years) that exist at the very earliest stages of star formation 
(Cabrit et al. 1990, Edwards, Ray and Mundt 1993). They occur not only in the massively active 
FU-Ori stages of star formation, but they also are to be found in the more quiescent, Classical T-Tauri 
Star (CTTS) stages. They appear to be generated by the interaction between protostars and accretion 
disks, within 10 stellar radii (R*) of the protostar (Hartmann 1992). 

What is the connection between bipolar outflows and protostellar jets? There is a growing consensus 
that bipolar outflows are a byproduct of protostellar jets (Snell et al. 1980, Shu et al. 1993), where 
the protostellar jet sweeps up ambient molecular cloud material into two thin shells, which manifest 
themselves as the observed bipolar lobes of CO emission. Once the molecular cloud material has been 
swept away (on a timescale of 10 5 years), the bipolar outflow disappears, leaving the protostellar jet to 
erratically fire away for a further 10 6 - 10 7 years. 

Our interest in protostellar jets is sparked by the following mass processing argument. Suppose 
the solar nebula formed a protosolar jet. Such a jet would have ejected 10~ 5 — 10~ 3 M of "rocky" 
material (i.e. all elements excluding H and He) from the solar nebula. If only 10% of this material 
were to fall back to the solar nebula then we would have 10~ 6 — 10~ 4 M of rocky ( possibly refractory) 
material being contributed to the solar nebula over a 10 6 — 10 7 year period. Given that the "rock" mass 
of the planets is of order 10~ 4 M , a protosolar jet may have made a significant contribution to the 
chemical structure of the solar nebula. Indeed, it is possible that ejecta from the jet flow may have been 
incorporated into the best preserved samples of the solar nebula: the chondritic meteorites. 

Of course, such an argument does not prove that protostellar jets formed chondrules, but only 
provides a plausible basis for constructing a theory of chondrule formation. To create such a theory, our 
first task is to investigate the thermal environment of the jet formation region. This is a function of the 
formation distance of a protosolar jet from the protosun. A distance which is uncertain, but model fits 
to the Li 6707 A, Fe 1 4957 A, Fe II 5018 A lines in FU Ori, plus the lack of extinction and infrared 
excess in the wind, suggests a formation distance of < 10i?* 0.1 AU for the solar nebula), where 
.R* is the radius of the protostar (Hartmann 1992). The temperature of the disk surface (T s ) at such 
distances, is given (approximately) by the formula (Frank, King and Raine 1992) 



' (M/M ) (M/1O- 7 M yr- 1 ) \ "\ 

A— —3 L (1 - (i?*Af /2 ) 7 



T s (r) = 850 K | ^ ^ >- | (l - (Rjr)^) ' , (1) 

where M is the mass of the protostar, M is the mass accretion rate, and r is the radial distance 
from the protostar's centre. The maximum temperature of the midplane of the accretion disk (T m ) can 
be approximately determined from the formula 

T m » T^r) 1 / 4 , (2) 

where rj is some number of order 1, and r is the optical depth (ibid., and Cassen 1993). Since r may 
have values as high as 10 4 , it is probable that at distances of order 0.05 AU from the protostar, the 
temperature at the midplane maybe high enough to melt or vaporise any rocky bodies that happen to 
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be in that section of the accretion disk. This raises the possibility that the protostellar jet winds may 
ablate droplets from the surface of rather warm "rock" bodies and if the droplets are small enough, they 
may be ejected from the inner accretion disk by the drag force of the jet flow. 

It may be asked, however, why embrace the idea that chondrules are ablation droplets? We are en- 
couraged to advance such a hypothesis, because chondrules appear to have been extruded or drawn from 
one or more extended magma bodies (Dodd and Teleky 1967). Also, an ablative process readily forms 
chondrule-like spheres from meteors (Brownlee et al. 1983). These "ablation-spheres" are chondrule- 
like in size and shape. They also share some physical similarities, e.g., they often contain relict grains. 
Despite this supportive evidence, there are, however, some significant problems to overcome before one 
can believe the ablation hypothesis. 

First, we require the appropriately sized rocky-bodies to enter into the jet formation region. Pro- 
tostellar systems generally have accretion disks and large bodies could be caught up in the inward 
accretional flow. However, the existence of chondrules in single ordinary chondrites that have (admit- 
tedly uncertain) 129 Xe*/ 127 Xe ages ranging over 10 7 years (Swindle et al. 1991, Swindle and Podosek 
1988), suggests that the accretional flow in the solar nebula was small enough to allow the survival of 
the chondrite precursors. Indeed, it is possible that such accretional inflow may have been dependent 
on height within the accretion disk, where only the gas at higher altitudes took part in the accretional 
flow (see §5). 

Another mechanism for radial migration is that of gas drag due to the velocity difference between 
nebula gas and the near Keplerian motion of large bodies around the protosun. The infall velocity for 
material subject to this gas drag is ( Whipple 1973, Weidenschilling 1977) 

dr ^ 2rAV 

dt VKeptD 

where Vx e p is the Keplerian velocity, AV is (~ lO^^Vxep) the difference between the gas angular velocity 
and Vxep, while to is the time scale for gas drag to influence the motion of the body. In the Stokes drag 
regime, 

where pb is the mass density of the body and r] g is the viscosity of the gas, and rj od is the radius of 
the "chondrule-f odder" body. Using the appropriate formula for the gas viscosity (e.g. Eq. (2.3.5) of 
Liffman 1992), one can obtain, from Eq. (3), the timescale for orbital decay ijdecay) 

r„ ~ lO 5 ( (W10 3 cm) Mlgcmj \ 

Tdecay ~ 1U ^ ((T^/lOO)^ K ) J ^ ' ^ 

with T g being the temperature of the gas. 

So, objects with a radius of around 10 m would fall into the Sun on a timescale of about 10 5 years. 
This analysis, however, ignores the change in the size of the body due to accretion of dust from the solar 
nebula. If the dust to gas mass ratio is of ~ 1, accretion of material will increase the size of the body 
and eventually stop its infall into the Sun (Weidenschilling 1988). On the other hand fragmentation 
due to inter-body collisions would have been a source of smaller material which may have eventually 
reached the boundary layer between the protosun and the solar nebula. 
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There may have been other mechanisms that brought macroscopic material into the inner regions 
of the solar nebula. We will simply note that it is a plausible assumption and one we require for our 
model. 

Once we have these large bodies in the jet-formation region, we are faced with the following scenario: 
the high speed jet-flow will probably occur at z > the scale height of the accretion disk, H, since the 
outflow jet must be governed by the conservation of mass equation, and so high velocities will only 
occur when p g (z) <C P#(0). On the other hand, our "chondrule-fodder" bodies will be located on or near 
the midplane of the nebula (z — 0). This separation between the critical points in the wind and the 
chondrule-fodder bodies would appear to be the death-knell for our ablation hypothesis - for how else 
can one produce ablation droplets if the hypothetical windflow is nowhere near our chondrule-producing 
planetesimals? 

To answer this, we examine protostellar jet theory. 

Most protostellar outflow models assume that magnetic fields provide the main coupling mechanism 
between the accretion disk and the outflow ( for recent reviews, see Bicknell 1992, Konigl and Ruden 
1993 plus Shu et al. 1993). In the most recent models, the magnetic driving mechanism involves the 
interaction between the dipole field of the protosun and the inner accretion disk (e.g. see Lovelace et. 
al. 1991, and Shu et. al. 1994). Studies of this protosun-nebula interaction can be traced back to the 
work of Freeman (1977), who assumed that the dipole field of the protosun was able to thread the inner 
(partially ionized) disk as shown in Fig. 1. In effect, the magnetic field of the protosun is "tied" to its 
surrounding accretion disk. 

In Freeman's model, the dipole field of the protosun rotates (approximately) with a rigid body 
velocity, so there exists a point away from the protosun (~ 0.04 AU), where the speed of the field 
sweeping over the disk equals the (approximately) Keplerian velocity of the disk. In Fig. 1, this position 
is denoted as the "synchronous orbit". For distances r greater than this synchronous distance, the 
magnetic field of the protosun has a greater angular velocity than the the accretion disk and so the 
magnetic field becomes "wrapped up" , i. e, the purely poloidal magnetic field of the protosun is converted 
into a toroidal field in the disk. 

In Fig. 1 , we show a side view of the resulting field structure. The toroidal field becomes the 
dominant field in the disk, and its direction reverses as one traverses the central plane. For such a 
configuration, the central plane of the nebula can become a "current sheet" where magnetic fields 
reconnect and the magnetic energy so released is converted into particle energy. 

As discussed in Priest (1994) the particle flow velocity (vf) obtained from the merging of magnetic 
fields (B) is close to the characteristic speed of a magnetic medium (i.e., the Alfven speed Ca), which 
in mathematical notation has the form 



B ( B \ ( 10~ 8 g cm 



3\ 1/2 

Vf « C A = = 3 ( — ^— ) ( — — km s" 1 , (6) 

where p g is the mass density of the gas and po is the permeability of free space. In the accretion disk, the 
direction of the reconnection flow would be roughly parallel to the central plane, with a small component 
in the z direction. 

To obtain ablative behaviour from gas flows with gas densities in the range of 10~ 6 to 10~ 8 g cm -3 , 
the required gas speeds range from 1 to 20 km s _1 (Liffman 1992). Comparing this to Eq. (6) suggests 
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that these "reconnection" flows may have the required densities and speeds to ablate the planetesimals. 
Some of the ablation droplets so produced would then, presumably, be caught up in the flow field that 
eventually becomes the main protostellar jet flow. 

It may turn out, however, that outflows are not powered by "wrapped- up" toroidal fields. Our 
purpose in describing this particular model is to demonstrate that, in the outflow region, the central 
plane of the accretion disk may be a highly active region, where conditions are conducive to chondrule 
formation by ablation. 



2. CHONDRULE COOLING 

A long standing problem in chondrule formation is how to explain the slow rate of chondrule cooling. 
Experimental simulation of chondrule formation suggests that chondrules cooled at a rate of 5 to 2000 
°C/hour (Hewins 1988). Such cooling rates are 3 to 5 orders of magnitude smaller than those expected 
for an isolated black body radiating heat directly into space. It would appear that chondrules were 
formed in a hot optically-thick medium or were produced in close proximity to other chondrules so that 
mutual radiation could damp the cooling rate. 

In our model, chondrules are produced in the hot, optically-thick midplane of the inner accretion 
disk around a protostar. They are produced by the ablative interaction between the initial stage of a 
protostellar jet wind and rocky chondrule "fodder" bodies that happen to stray into the jet formation 
region. If the subsequent droplets are small enough, they will be swept up by the wind and begin to 
move with the jet flow. 

As the jet flow moves away from the midplane of the accretion disk, it is likely that the gas in the 
flow will cool, since work is being done to expand the gas. Consequently, if the molten droplets are 
close to thermodynamic equilibrium with the surrounding gas flow, they must also cool at just about 
the same rate as the expanding gas flow. 

Naturally, to model this process we require a model of a protostellar jet. Unfortunately, a com- 
prehensive theory of protostellar jets is unavailable at this time. To partially circumvent this difficulty, 
we adopt a simple parameterised model of an outflow. First, we describe our system with cylindrical 
coordinates, the plane z = being the midplane of the accretion disk with the protostar residing at 
the origin (see Fig. 11). With this coordinate system, the steady state form of the continuity or mass 
conservation equation (V • (p 9 v ff ) = 0) is easily solved, for an axisymmetric flow, to give the equation 

2p g (z)Vg Z (z)irr(z) 2 = constant = M Q , (7) 

where p g (z) is the "r-averaged" gas mass-density as a function of z, r(z) is the cylindrical radius 
of the outflow, Vg Z (z) is the z component of the "r-averaged" gas velocity, and M is the mass- loss rate 
of the outflow. The factor of two in Eq. (7) arises, because the protostellar jet is produced from both 
sides of the accretion disk. 

To model the gas density and gas velocity profile within the disk, we assume that the disk constrains 
the radial size of the protostellar jet. If we now let 

/ Z\ \ m 

Pg(z) « Pg{zi) (— I , Z < Z X , 771 > , (8) 
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where z\ is the distance above the disk midplane where we can best define the flow variables. We can 
substitute Eq.(8) into Eq.(7) and obtain 



v gz (z) w v gz {z x ) » z < Zi , m > . 



(9) 



These solutions break down as z — > 0, this being the price we pay for not solving the momentum and 
energy equations. Equations (8) and (9) are only a first order fit to the mass conservation equations. 
However, they do, approximately, satisfy the density-velocity relationship observed in thermally driven 
winds, i.e., as p g decreases, v gz increases. 

If the gas pressure (P g ) can be globally modelled as a polytropic gas, then 

P g = Kg], , (10) 

where k and 7 are constants (7 is the ratio of the specific heats if the jet is adiabatic). Using the ideal 
gas law with Eqns (8) and (8), we can obtain an expression for the gas (and chondrule) temperature 
(T g ) as a function of z. 

/7, \ mty— 1) 

T g {z)^T g { Zl )[j) . (11) 

The gas flow and entrained chondrules decrease in temperature as they travel away from the central 
plane of the accretion disk. 

If we assume that the particles are entrained with the flow, and they start from z — 0, then we can 
use Eqns (9) and (11) to give the gas and chondrule temperature as a function of time (t). Solving for 
t in Eq. 9 gives 

1 l/(l-m) 

0<m<l, (12) 



z 



r-i 



where T\ = zi/((l — m)v gz (zi)). Substituting the above equation into Eq.(ll) gives, 



r-i 



m("f— 1) 



TM = TM)[j\ W 0<m<l. (13) 

The gas and chondrule cooling timescale is determined by the value of T\ and T g (zi). As we shall show 
in the next section, for T g (zi) ~ 1200 K, then z\ ~ 10 4 — 10 5 km and v gz (zi)) ~ 0.1 — 1 km s _1 . As a 
consequence, a minimum parameterization for T\ is 

= in 4 q (^i/10 4 km) 

1 (l-m)(^ 2 (z 1 )/lkms- 1 )- 1 ) 



The general decrease in gas temperature, and chondrule temperature, as a function of z within or 
near the accretion disk has important consequences for the collisional interactions between chondrules. 
Indeed, a simple analysis of this phenomenon leads to a model which can possibly explain the observed 
structure of compound chondrules and leads to a prediction for the physical structure of triple compound 
chondrules. 
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3. CHONDRULE COLLISIONS 
3.1. Hover particles 

Let us consider a chondrule that has just been ablated from its parent body. We assume, that 
the chondrule will be accelerated by the gas jet so that its motion is, initially, in the z direction. The 
equation of motion for the particle, parallel to the z axis, is given by 

.. C D 22 GMm p z 

m pZ = ~yPMz - z) na p - - ^ 2)3/2 , (15) 

where Cd is the coefficient of gas drag, while m p and a p are the mass and radius of the particle, 
respectively. 

Suppose the particle reaches a state where the gas drag is balanced by gravity, so that z — 0, and 
z = 0. Our equation of motion becomes (for z <r) 

n~ Cr > 2 2 GMmz h 

~ ~YP 9 V a P ^3— > ( 16 ) 

where Zh is the value of z for the hovering particle. For z close to the midplane, we should expect that 
v gz will be less than the sound speed, so if the mean free path of the gas (I) satisfies the relation 



l/2a p >10 (17) 



then Cd has the Epstein (1924) form 



C D « " ^ , (18) 
3v gz y nrrig 

where k is the Boltzmann constant and m g is the mass of a gas particle (in this case we assume 
monatomic hydrogen). Using Eqns (16), (18) and (7) we find 

/8fcT g \ 1/2 Mr 
Zh ~ V ) 2nGMp p a p ^ ' 

or 

(M/M )(p p /1 gm cm- 3 )(a p /0.1 cm) V ; 

Thus, Zh ~ 10 -4 AU <§C 0.1 AU ~ r as is required from the assumption leading to Eq. (16). Now, Zh is 
inversely proportional to the radius of the particle (a p ), so smaller particles will have higher hover heights 
than larger particles of the same density. To reach their respective hover altitudes, the particles will 
move relative to each other and may undergo collision interaction. If the collisional velocities are small 
enough and the temperatures high enough, we may have chondrules fusing together to form compound 
chondrules. 
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3.2. 



Compound Chondrule Formation 



Studies of compound chondrules have classified them into two general types: enveloping , where 
one chondrule envelopes the other and adhering , where one chondrule forms a "bump" on the other 
(Wasson 1993). In this paper, we restrict our attention to adhering compound chondrules. 

Most adhering compound chondrules consist of a smaller chondrule stuck to a larger particle, where 
the smaller chondrule was plastic at the time of collision (Wasson et al. 1995). This fact immediately 
poses a major problem for most chondrule formation theories which use simple "flash" heating scenarios. 
Since, if the chondrules were all formed in the one flash heating event, then we should expect the larger 
chondrules to be plastic at the time of collision. 

In the Jet model, chondrules are continuously formed by an ablative process as chondrule fodder 
bodies stray into the jet formation region. As discussed in the derivation of Eq. (20), it is reasonable to 
assume that 0.1 cm particles will be supported by the jet flow at a z distance of 10~ 4 to 10~ 3 AU above 
the midplane of the accretion disk. Let us suppose that at such distances, the temperature of the gas 
in the Jet stream is around 1200 K - the approximate solidus temperature of a chondrule. 

Assuming a fairly steady jet flow, a hover particle in such a position should have sufficient time 
to equilibrate with the gas temperature and be relatively non-plastic. The hover particle will only see 
smaller particles flying past it in the jet flow, since larger particles of approximately the same mass 
density will have a lower hover altitude. 

As these smaller particles fly past the hover particle, the thermal inertia of these particles will give 
them higher temperatures, and therefore greater plasticity than the hover particle. If the hover particle 
collides and fuses with these particles, then the observed structure of adhering compound chondrules 
will be obtained (for a schematic depiction of this process, see Fig. 2 ). 

To give this idea a quantitative context, we note that the equation for the rate of temperature 
change of a spherical chondrule has the form 



(Liftman 1992 and references therein). 

The left hand side of the equation describes the time rate of change of the heat energy of the body, 
where a p is the radius of the particle, p p is the density of the particle, Cy is the specific heat per mass 
of the body (~10 7 erg g _1 K _1 for chondrules), Tg is the surface temperature of the body, and t the 
time. The first term on the right hand side describes the energy added to the body by the gas/body 
interaction, where A is the heat transfer coefficient and q is the gas/body heat transfer rate per unit 
surface area. The last two terms are radiation terms, where T e is the radiation temperature of the 
surrounding environment, e a / e is the absorptivity/emissivity of the surface (absorptivity and emissivity 
are assumed to be the same in this case), and a is the Stefan-Boltzmann constant. In the following 
discussion, we will assume local thermodynamic equilibrium applies and that T e = T g . 

For chondrule formation, the size of the chondrules (< 1 cm) and the expected low gas density 
(< 10~ 6 g/cc) indicate that the mean-free path of the gas is large relative to the size of the chondrules. 
For this "free molecular" flow regime Awl and q has the form 




dT s 



dt 



,q + A-Ka 2 p e a aT^ - 4ira 2 p e e aTg 



(21) 



q = Pg\Vg - V p \{T r 



rec 



T S )C H , 



(22) 
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SIDE VIEW 








Fig. 1. — The dipole field of the protosun interacts with the accretion disk, such that for distances greater 
than a critical distance, (~ 0.04 AU) from the protosun, the poloidal field of the Sun is "wrapped-up" 
by the accretion disk. This process may cause the central plane of the inner accretion disk to become 
a "current sheet" produced by the merging magnetic field lines. This merging will be due to magnetic 
diffusion and the switch in the sign of the toroidal field as one traverses the central plane. 



< 1200 K 



1200 K 



> 1200 K 




rising plastic chondrule 



Fig. 2. — Schematic depiction of compound chondrule formation in a jet flow. The gas flow is depicted 
within the accretion disk, where the gas temperature decreases as z (the distance from the midplane of 
the disk) increases. A "large" hovering chondrule, with a temperature near or at the solidus temperature 
of 1200 K, collides with smaller chondrules that are entrained in the gas flow. Due to thermal inertia, the 
smaller chondrules are warmer and more plastic than the hover chondrule. In this way, one can obtain 
the observed phenomena that smaller chondrules were plastic at the time the compound chondrule was 
formed. 
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(Hayes and Probstein 1959, Probstein 1968) where p g is the mass density of the gas, v g is the gas velocity, 
v p is the particle velocity, T rec is the recovery temperature, and Ch is the heat transfer function for free 
molecular flow. T rec and Ch have the forms 



T 

VPC 



7 + 1 

and 



2 7 — I 

2 7 + 2( 7 - l)s gp - a5 + a _ 2 + S9p7r -o.5 exp( _ s 2 p)/erf(Ssp) 



(23) 



Ch = -^[^s^i-sl) + (0.5 + si) eri(s gp )i (24) 

(i&id. ) where T g is the gas temperature, m g the mean gas particle mass (~ 1.66 x 10~ 24 g), 7 is the 
ratio of specific heats, k is Boltzmann's constant, erf(s) is the error function (= ^= exp(—t 2 )dt ), and 
s gp is the ratio of the relative streaming gas speed and the most probable Maxwellian gas speed, 



\v g ~ v p\ 



m 



(25) 
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We are considering flows at the very base of the protostellar jet, so we should expect that s 9J) < 1. 
Thus, q becomes 

,»*(T.-T.)(I±l)^/g. (26) 



9 



or 



q » 5 x 10 7 ( p -) ( ( -^-) - ( ( J ergs" 1 cm" 2 s" 1 . (27) 

H V 10 s gem -3 / VVlO 3 K/ V10 3 KJJ V7-I/ V 10 3 K & v ; 

When p 9 >, 10 -7 g cm -3 , 7 = 5/3 (monatomic gas), and T g ss 1200 K, then the q term in Eq. (21) 
dominates the radiative terms, since 

(tT 4 « 5.7 x 10 7 f-^-V ergs" 1 cm" 2 s -1 . (28) 
y V10 3 K/ 

So Eq.(21) becomes a first-order linear differential equation in T s , and has the solution 

T s (t) » T fl + (T s (0) - T s )exp(-t/r) , (29) 

where r has the form 



a p PpC v /7 - 1 



37rm 3 



-^ItTTiVwf' (30) 



or 

^ (a p /0.1 cm)(p p /3 g cm- 3 )(CV/10 7 erg g" 1 K" 1 ) / T s x 
T ^ 5 (p fl /10"8 gem" 3 ) \Wk) S [6L) 

So, if we have p g >, 10~ 7 g cm -3 then particles with a radius between 0.01 cm and 0.1 cm will have 
a minimum temperature equilibration timescale of 0.5 to 5 s. Similar timescales are obtained for the 
low density case (p g <C 10~ 7 g cm -3 ), where the radiative terms in Eq.(21) dominate the "g" convective 
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term. In this case, Eq.(21) can only be solved numerically. The results from this calculation are shown 
in Fig. 3. 

As can be seen, The e-folding timescales in the radiative case turn out to be about the same as 
those for the pure conduction case. One should note that these are minimum timescales, because we 
have neglected the latent heat of fusion and the diffusion time for heat to travel from the surface of the 
particle to its centre. 

Despite these caveats, it is probable that the timescales for thermal equilibration are not much 
greater than as stated in Eq.(31). In such circumstances, we can only expect the smaller chondrule to 
be plastic at the time of collision if, in the collision zone, there is a sharp decline in gas temperature. 
We suggest that a temperature drop of 10 to 100 K s _1 is required for compound chondrules to form. 
We also expect that the collision zone for compound chondrule formation must be quite thin, since the 
maximum time that the smaller secondary can remain plastic is ~ 20 s and the maximum collision 
velocity is ~ 0.1 km s -1 ( a higher collision velocity would fragment the chondrules, see Vedder and 
Gault (1974)). These numbers give a maximum length scale for the thickness of the compound chondrule 
formation region of around 1 km. 

Now, suppose we have a chondrule hover-particle in a region of the flow where the ambient temper- 
ature of the gas flow is much less than the solidus temperature of the particle. If this solid hover-particle 
collides with smaller particles which are entrained in the jet flow, and if the relative speed of the two 
particles is greater than 0.1 km s -1 then it is likely that one or both of the particles will undergo some 
damage from the collision. Depending on the relative speed of collision, this damage can range from 
slight chipping to complete fragmentation. 

It is well known that many chondrules have undergone chipping and fragmentation. In the Jet 
model, one can provide a cause for this damage: collisions between chondrules. One can also give a site 
for where this damage will take place: in the jet flow at higher (and hence colder) altitudes than the 
compound chondrule formation region (see Fig. 11). In this way, the jet model can provide a unified 
scheme linking both compound chondrule formation and chondrule fragmentation. 

Of course, like all theoretical mutterings in this field, the above results and ideas should be (and 
will be) treated with a reserved caution. A physical theory is of little use, unless it can be used to 
predict as well as to explain. So far we have attejrn^ed^the latter, in the next section we will try the 
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Fig. 3. — The radiative cooling of 3 particles of different radii from 1250 K to 1200 K. 
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former. Let us throw caution to the winds (no pun intended) and predict the physical structure of triple 
compound chondrules. 



3.3. The Weather Vane Effect 

If compound chondrules were formed in protostellar jets, then gas drag would orient the compound 
chondrules immersed in the jet flow. This is simply because the spherical symmetry of the simple 
chondrule is destroyed. A secondary chondrule will act like the tail on a weather vane, orienting a 
binary compound chondrule such that the smaller secondary will be on the "downstream" side of the 
primary. A subsequent collision and fusion with another chondrule will produce a triple compound 
chondrule where the secondary chondrules are separated by a minimum avoidance angle. This situation 
is shown schematically in Fig 4. 

ttt 

Jet Flow 




incoming secondary 



Fig. 4. — The "Weathervane" Effect. Secondary chondrules in a streaming gas flow, will tend to orient 
the compound chondrule such that the secondary is pointing in the direction of the flow. Any incoming 
chondrule will hit the rear end of the compound chondrule. Hence the prediction that secondary 
chondrules in a triple compound chondrule will tend to avoid each other. 

This simple scenario does have some initial complications, the major being that if the impact 
parameter (the offset distance between the centres of the particles, see Fig 4) between the colliding 
chondrules is greater than zero, then the binary compound chondrule will be set spinning. This rotational 
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motion will, however, be damped, because the secondary chondrule will be travelling half the time with 
the gas jet, and half the time against the jet. If the secondary chondrule has a cross-sectional area A, 
then the rate of change of the rotation rate, u, is given by 

Cd A Cd A 
ICo ~ —p gaps ( Vgz - a psUJ f— —p gClps ( Vgz + a ps u) 2 — , (32) 

where / is the moment of inertia of the system, and a ps is the distance between the primary and 
secondary chondrules. If we substitute the Epstein solution for Cd (see Eq. 18) , we obtain 



^imv*f±) Ut (33 ) 

3 y nrrig \ I / 

which has the solution 

lu w u exp(-t/r sd ) , (34) 
where ujq is the initial rate of rotation, and r s d - the e-folding "spin-down time"- has the form 

Tsd ~ 8p 9 Aal s y 8kT 9 ■ [6b) 
To obtain an estimate for t S( i, we note that 

I ~ (m p + m s )a 2 ps < 2m p a 2 ps , (36) 

and 

A ~ 7ra 2 s , (37) 

where m p (a p ) and m s (a s ) are the masses (radii) of the primary and secondary chondrules, respectively. 
Combining Eqns (36) and (37) with the observation that a s >, 0.1a p and putting it all into Eq.(35) gives 



(p g /l0- 8 g cm- 3 ) W T< 



9 



So for a range in gas densities of 10 6 to 10 10 g cm 3 the spin-down timescales for a binary compound 
chondrule range from around a minute to about a week. 

Once a binary compound chondrule has stopped spinning, then it will take on a specific orientation, 
where the smaller secondary will be on the downstream side of the primary chondrule. As shown in 
Fig 5, the system can still oscillate between two extreme positions, but there will always be a section 
of the primary chondrule which will be shielded from colliding with particles that are entrained in the 
gas flow. 

From Fig 5, one can deduce a minimum avoidance angle 6 that two secondary chondrules on a 
triple compound chondrule should obey if the initial binary compound chondrule has stopped spinning. 
Given that the distance between the centres of the primary and secondary chondrules is a ps = a p + fa s , 
where / is the deformation factor of the secondary chondrule, then 

cos (fl) = a p-* s = 1 -^/a P (3g) 
a p + fa s l + fa s /a p 
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A plot of a solution to this equation is shown in Fig 6. If the weather vane effect did not occur 
then we would expect a uniform distribution of angles between secondaries for each value of a s /a p . If 
the weather vane effect does occur, then the distribution should be skewed towards the 'allowed angles' 
shown in Fig 6. This is what we see in the few data we have been able to scour from the literature. 

Again, we should make clear that the weather vane effect will only be observed if the timescale for a 
binary compound chondrule to stop spinning is less than the chondrule-chondrule collisional timescales. 
It will also only be observed in triple compound chondrules (i.e. one primary and two adhering sec- 
ondaries). Compound chondrules that have two or more adhering secondaries will have unpredictable 
orientations and so an additional secondary may land anywhere on such a compound chondrule. 

Finally, we should point out that we have ignored gas flow fluctuations that are parallel to the 
central plane of the accretion disk and that the derivation of Eq.(33) assumes v gz remains constant 
during one revolution of the spinning compound chondrule. This latter assumption is probably true 
just after the formation of the compound chondrule, but may not be true as the compound chondrule 
stops spinning. 

To see this, we note that a compound chondrule will stop spinning when it receives enough torque 
from the gas jet to remove all the angular momentum and reverse the chondrule's initial spinning 
direction. The length of time the secondary will be moving against the gas jet is n/ou. The torque is 
approximately Cr>p g v 2 gz Aa v l1. So the minimum rotation rate cu m i n , is given by 

C 7T 

IuJmin ^ —p g v 2 gz Aa p . (40) 

Thus 

^ ^ CD "f Aa * 'f. (41) 

Using the Epstein drag law (Eq.(18)) plus the approximate values for / and A (Eqs (36) and (37), 
respectively). One can show that 

^ „ . 015 ^^ ( (^/°.lk.ns-')(V10-n gcm - 3) ) V, ^ ^ 
{a p /0.1 cm.y V (Pp/3 gem 6 ) J 

where we have assumed that T g = 1000 K and that m g = m H . So the period of rotation (2tt/u) of this 
slowly spinning compound chondrule is approximately 400 seconds. In a 0.1 km s -1 gas flow, this gives 
a scale length (L) of order 10 to 100 km. 

So the relevant Reynolds number (Re) of the flow (see Liffman 1992) is given by 

Re « 8.7 x 10 3 (L/10 3 km)(^/0.1 km s* 1 )(p 9 /10~ 11 gcnr 3 )^/(T 9 /10 3 K)(m H /m g ), (43) 

so, with Re ranging from 100 to 1,000, v gz may undergo fluctuations which may modify the spin-down 
time of Eq. (35). 



3.4. Chondrule Reheating 

If we combine the theory that we have acquired from §2 & §3, we can present a scenario for the 
phenomenon of chondrule reheating, where high temperature rims are observed around chondrule cores 
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Fig. 5. — The extreme positions of the secondary chondrule while oscillating from side to side. 
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Fig. 6. — Comparison of theory vs observations for triple compound chondrules. We plot the angle 
between the centres of secondary chondrules vs the radius ratio of the largest secondary to the primary. 
The largest secondaries were used for a s values, as they would show 'disallowed' chondrules, if any were 
present. As can be seen, all the observed angles lie in the 'allowed' zone. The line separating the 
'forbidden' and 'allowed' zones is given by Eq.(39). 
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(Kring 1991). We suggest that reheating may be a consequence of a chondrule overshooting its hover 
altitude and then oscillating, in a damped manner, around the hover altitude. Each time the chondrule 
decreases its altitude, it will undergo an increase in temperature. If the altitude decrease is large enough, 
the corresponding increase in temperature will be sufficient to remelt the chondrule. 

The actual mechanism whereby these high temperature rims are formed will not be discussed in 
any depth here. We only suggest that as the chondrule oscillates around its hover point, it accretes 
dust-grains/molten droplets and that this material forms the foundation of the rim once the chondrule 
is reheated. 

To understand how these oscillations arise, we note that (see also §3.1) the equation of motion for 
a particle, parallel to the z axis, is given by 



m n z 



C D 2 2 GMrripZ 

-7rP 9 {v gz - z) na p - , 2 , 2W2 • (44) 



> - ) r 'l" ' ' ; ""/' ( z 2 _|_ r 2)3/2 

As noted in §3.1, the drag coefficient (Cd) is likely to have the Epstein form 



C D « " /?*?k. (45) 



3v gz V irrrig 

Combining Eqs (44) and (45) and assuming z C r, we have: 

where z h is the hover height given by Eq.(19). 
Let x = z — z h , then Eq.(46) becomes 

rripX + cx + kx = 0, (47) 

where 

/ 8nkT g \V* 

o - Ta pPg 

and 



GMm D , . 

k = — (49) 

If we assume that c and k are approximately constant then c 2 < 4km p gives solutions of Eq.(49) that 
oscillate, with decreasing amplitude, as a function of time (t). If c 2 = Akm p or c 2 > Akm p then the 
solutions of Eq. (49) are critically or strongly damped and the particle approaches its hover height 
without oscillation. 

The case c 2 = 4km p translates into an equation which gives the critical density for particle oscillation 
(Pgc)- 

12 (op/0.1 cm)(p p /l gm cm*){M/M & yl\m g /m H yl* 3 
P 9 c-2.7xW (r/0.1 AU)3/2(T/1Q3 gCm ' (5 ° j 
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If p g < p gc then c 2 < 4km p and the particle will oscillate around the hover point Zh- It is our contention 
that it is this oscillation process that causes chondrule reheating. Since the temperature of the jet flow 
increases as z decreases, thus as the particle travels to lower altitudes its temperature must increase. 

To be specific, when p g < p gc then Eq. (47) admits the solution 

z = ^ + fM e - w sin/it, (51) 
p 

where z(zh) is the z speed of the chondrule when it first reaches z = Zh, 

A* 2 — ~— 7~^r— > (52) 



and 



Thus the period of the oscillation is 



4m 2 



2m p 



(53) 



H ~ (Mo) 1 / 2 
while the damping timescale is given by 



(r t\T , days (54) 



1 (flp/0.1 cm)(p p /l gm cm 3 )(m g /m H ) 1/2 , , , 

b~ (p 9 /10- 12 gmcm-3)(T 9 /10 3 K)V2 ' days ^ 

and the amplitude of the oscillation has the form 

iM = 1 . 6xl0 *(£M/^ km ,56) 

Comparing Eqs (56) and (20) shows that if z[zh) is comparable to the gas flow speed, then the amplitude 
of the oscillations can be a significant fraction of the hover height. The corresponding temperature 
variations may also be significant. 

Of course, from Eq.(48), c is a function of p g and T g and cannot be considered a constant, so 
the above analysis should only be treated as an informative approximation to the complete system. A 
slightly more realistic analysis can be obtained by constructing a computer simulation of the flow. We 
were able to accomplish this by assuming a z velocity for the outflow wind of the form 



v g {z) - 



/ z \ m 

0.1 - kins" 1 , z < z h , (57) 



with m = 0.1, with the orbital (angular) velocity of the gas assumed to be Keplerian. The density 
structure of the gas flow could then be easily deduced from Eqs (7) and (57), and had the form 



M/lO" 8 M /yr 

pJz) w 4.5 x 10' 12 / — ) 1 °/ ; ; g cm" 3 , 58 

l gy! (v g (z )/0.1 km s _1 r/0.1 AU 2 s ' y ; 
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while the temperature of the gas was given by 

T g (z) « 1200 (^-jj K, (59) 

where we set 7 = 5/3. 

As can be sen from the above formula, T g (z) will be greater than 1900 K when z £ 0.001^. In 
such cases, we simply set T g (z) = 1900 K. Finally, to compute the temperature of the particle, we did 
a full time integration of Eq.(21). 

In Fig. 7a we show the motion of a chondrule-like particle (a p = 0.1 cm, p p = 3.8 g cm" 3 ) released 
from very close to the midplane (z = 100 cm) into a jet flow with a total mass loss rate of 10~ 9 M yr _1 . 

The particle was released at a distance (r) of 0.05 AU from the center of a solar mass protostar, 
where 0.05 AU is also the assumed radius of the protostellar jet. The particle underwent a series of 
damped oscillations the period of which was slightly less than 6 days, while the damping timescale 
was approximately 30 days. The corresponding temperature of the particle is shown in Fig. 7b. The 
particle temperature has a damped "saw-tooth" pattern, the maximum temperature of the resulting 
temperature peaks being approximately 1420 K. Of course higher temperature values can be obtained 
if the temperature gradient is steeper than what we have assumed (e.g., if m > 0.1). 

While this oscillating motion is of interest, the particles can also undergo a critically damped 
trajectory. In particular, if p g > p gc then Eq. (47) gives the solution 

z = Zh + J^L( e (-H-V=^)* _ eM-v^)*), (60) 
2y/-fi 2 

The trajectory and temperature of such a particle are shown in Figs 8a & b, where now the mass loss 
rate of jet is assumed to be 2 x 10~ 8 M yr -1 . 

From this simple model, we would suggest that chondrule reheating occurs when the Jet flow is 
declining in mass flux, i.e. , during the later stages of a CTTS's evolution. 



4. THE CHONDRULE SIZE LIMIT 

Why do chondrules have an upper size limit? Possible solutions such as size-limited precursor dust 
balls to aerodynamic sorting have been suggested. In the Jet model, one has to eject particles from the 
inner accretion disk and ensure that they arrive in the outer parts of the accretion disk. Clearly, there 
must be some size limit to this process and this size limit must be dependent on the "strength" of the 
jet flow, which in turn is dependent on the gravitational potential of the inner accretion disk. To derive 
a quantitative relationship between chondrule and accretion disk, we must first describe the behaviour 
of a droplet in a flow. 

Suppose we have a droplet of molten material that is subject to a streaming gas flow. The surface 
drag of the flow will tend to "rip" the droplet apart, while the surface tension of the melt will try to 
minimise the exposed surface area of the droplet and keep the droplet together. The balance between 
these two conflicting forces produces a stable droplet of maximum radius a p , the formula for which is 

C DpgVl p 



a p ~ 7< I 72 ' ( 61 ) 
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Z vs Time Temperature vs Time 
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o 




Fig. 7. — (a) A chondrule-like particle is released into a protostellar jet flow at a distance of 0.05 AU 
from a protostar. The mass loss rate of the flow is set at 10~ 9 M yr _1 . Such a low mass loss rate 
ensures the gas density is lower than the critical density given in Eq. (50). In such a circumstance the 
particle will oscillate, in a damped manner, around the hover point (see Eq. (51)). (b) The temperature 
of a particle as it undergoes the trajectory shown in Fig. 7a. 



Z vs Time Temperature vs Time 




Fig. 8. — (a) A chondrule-like particle is released into a protostellar jet flow at a distance of 0.05 AU 
from a protostar. The mass loss rate of the flow is set at 2 x 10 -8 M yr _1 . Such a low mass loss rate 
ensures the gas density is higher than the critical density given in Eq. (50). In such a circumstance 
the particle will approach the hover point in a critically damped manner (see Eq. (60)). (b) The 
temperature of a particle as it undergoes the trajectory shown in Fig 8a. 
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where 7 is the surface tension of the molten material, and Weo is a dimensionless factor called the 
"critical Weber number" . We accounts for the non-uniformity of the gas drag pressure over the surface 
of the droplet. Typically, We « 10 (Bronshten 1983). 



The streaming gas flow will not only determine the stable size of the droplet, but it will also subject 



the particle to a drag force given by ^-p g v 2 A p , where A p is the cross sectional area of the particle that 



is facing the gas flow. 

One can show (Liffman and Brown 1995) that the work done, W, by the wind - flow in ejecting a 
particle is 

W = ^fA p < Pg v% > L, (62) 

where C* D is a "representative" value of the drag coefficient during the propulsion phase, < p g v 2 gp > 
is the mean value of p g v gp , and L represents the actual length of the propulsion stage. 

For a particle to escape the disk, we require that 

2 kGMm„ 

< p ' v "> >= RCIA-L^ (63) 

where k is > 1, G is the gravitational constant, m p is the mass of the particle, and R is the initial radial 
(or semi-major axis) distance of the particle from the protostar. Substituting Eq.(63) into Eq.(61) gives 

We Q C* Dl RA p L 

ap ~ kC D {0)GMm p ' [b ] 

where Cd(0) is the value of the gas drag coefficient at the position where the particle is formed. If we 
assume that our particle is spherical (as are most unfragmented chondrules) and that our propulsion 
distance L is proportional to the height of the disk at a distance R away from the protostar, i.e. L = \R@, 
then 

-3AVye C|,7# 1+/ Y /2 



dp ~ 



4C D (0)kGMp p 



(65) 



Now Cd, generally, decreases with increasing gas-flow speed. So, we should expect that C* d /Cd{Q) ~ 
since C* D is an 'average' gas drag coefficient over the entire propulsion stage, while Cd(0) samples the 
gas flow at the beginning of the propulsion stage, where the gas flow is probably at its slowest. Setting 
A ~ 0.01, (5 ~ 1, (typical, approximate values for the disk height), k ~ 1,M M ,and We Q ~ 10, we 
obtain 

R (AU). (66) 



a p (cm) <, 0.4 



1 



The values for the surface tension, 7, and the density, p p , are material dependent. We are interested 
in Fe-Ni and silicate chondrules, so we shall use the surface tensions of meteoric iron and stone, which 
are 7i ron = 1,200 and 7 s t one = 360 gs~ 2 ( Allen et al. 1965) with corresponding mass densities of 
Piron = 7.8 and p s tone = 3.4 gcuT 3 . Substituting these values into Eq.(66) gives us the approximate 
radius of iron and silicate droplets as a function of distance from the protostar: 

r iron (cm) < 5R (AU) and r stone (cm) < AR (AU). (67) 
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If chondrules were formed from this ablative process, then to obtain the observed maximum chon- 
drule sizes we require that R <, 0.1 AU (and k > 1 ). This is consistent with the conclusion that 
protostellar jets are formed in an accretion disk at a distance of 0.05 to 0.1 AU from a solar-type 
protostar (Camenzind 1990, Hartmann 1992). 

Although this result is encouraging, there is a contradiction between the ideas presented in §3 and 
§4. In §3, we require the chondrules to remain stationary in the flow so to allow the formation of 
compound chondrules, while in this section we have shown that the chondrule size limit is determined 
by the "power" of the jet flow to eject particles from the inner accretion disk. This contradiction would 
appear to severely limit, perhaps destroy, our theory. After all, how can one require particles to be 
stationary relative to the accretion disk and also expect them to be ejected at speeds > the escape 
speed of the protostellar system? 

The only way out of this contradiction is for the jet flow to be highly variable in both density and 
velocity. Compound chondrules would presumably form when the flow is relatively quiescent, while 
chondrule ejection would perhaps occur when there is a major increase in the density and/or velocity of 
the flow. Observations do suggest that jet flows are highly variable in their behaviour (Edwards et al. 
1993). For example, Mundt (1984) obtained observational evidence of variations in young stellar winds 
that vary on time scales of months. Despite this observational support, the resolution of this problem 
must await a coherent theory of protostellar jet formation, which in turn must await high resolution 
observations of protostellar jets. 

Besides determining a maximum size limit for chondrules, the ejection of chondrules by a proto- 
stellar jet also, indirectly, causes chondrules and chondrule fragments to be size sorted. 



5. SIZE SORTING AND DISK ACCRETION 

It has been known for many years (Dodd 1976) that chondrules, both silicate and metal, are size- 
sorted in meteorites, i.e., particles in a particular chondrite satisfy the relation (Skinner and Leenhouts 
1993) 

Pgdp « constant . (68) 

It has recently become apparent that size-sorting also applies to fragments of chondrules (Skinner and 
Leenhouts 1991). This shows that chondrites act as "size-bins". Chondrules and their fragments were 
formed, mixed and later sorted by size within the chondrite forming regions of the solar nebula. 

The usual explanation for this phenomenon (e.g. Dodd 1976 and references therein), is some form of 
aerodynamic drag effect, where larger and/or denser particles can travel further into a resistive medium 
than can smaller and/or less dense particles. We too shall use this idea, for it arises naturally from the 
Jet model. 

To see this, suppose that a particle is ejected from the inner accretion disk at speeds close to 
or exceeding the escape velocity. Suppose further that the angular momentum of the particle is high 
enough to eject the particle from the gas flow and allow it to move across the face of the accretion disk. 
If the particle is not subject to gas drag, from the upper-atmosphere of the accretion disk, it will simply 
move out of the system and into interstellar space. If the particle is subject to gas drag then, given 
sufficient drag, the particle will be recaptured by the protostellar system and may fall into the outer 
parts of the accretion disk. 
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A stream of such particles will be size-sorted. Since the smaller, less dense particles will fall close 
to the protostar, while the larger, more dense particles will fall further away from the protostar. Thus, 
we should expect that meteorites from the outer parts of the solar system will contain larger particles 
than meteorites from the inner solar system, and indeed, chondrules from the carbonaceous chondrites 
are (usually) larger than chondrules from ordinary chondrites. 

There is an additional consequence of this model which has to do with angular momentum transfer. 
A particle in Keplerian orbit around a star has an angular momentum that is proportional to the square 
root of the distance between the particle and the star. Thus, a particle that is ejected from the inner 
accretion disk, and subsequently stopped by gas drag, will simply fall back into its original orbit unless 
angular momentum is transferred from the disk to the particle. If gas drag allows the disk to transfer 
angular momentum to the particle, then the particle can fall to the outer parts of the disk. 

Clearly, if the disk gas gives up angular momentum, then it must move in towards the protostar. 
Since protostellar jets are fueled by disk accretion, one obtains the schematic picture of a jet flinging out 
material into the disk such that the disk will accrete onto the protostar and refuel the jet. Heuristically, 
one can think of chondrules and other associated Jet ejecta as delayed protostellar Jet fuel. 

To turn this qualitative speculation into quantitative speculation, we have to model the particle 
ejection process. Unfortunately, this is a difficult thing to do, since there is no consensus on how 
protostellar jets work. So, we make the following assumptions: 

(1) We suppose that our chondrule particles are, initially, in a circular Keplerian orbit of radius 
R (< 0.1 AU) from the protostar. 

(2) We assume that the protostellar jet gives the particles an initial "boost" velocity i(0) that 
is comparable to the protostar 's escape velocity at that point. Our tentative justification for this 
assumption is that protostellar jets are observed to have speeds comparable to the maximum escape 
velocity of a protostellar system. So, we presume that particles initially entrained in such a flow may 
also obtain similar speeds. 

(3) Finally we assume that the particles have, initially, no radial velocity, and that the self-gravity 
of the disk is negligible compared to the gravity of the protostar. The former of these two assumptions 
comes from the observation that jet flows tend to be perpendicular to their respective accretion disks. 
Protostellar jets do have a nonzero radial velocity, but we have, for simplicity, ignored this component. 
The validity or otherwise of these ideas is discussed at some length in Liffman and Brown (1995). 

Given these assumptions, the equations of motion for a particle become 

,_*- (69) 



and 



v e = rO=-, (70) 
r 

.._ GMz 

where h is the specific angular momentum of the particle and has the value 

h = VGMR . (72) 
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The value of h is a constant, since it is assumed that there are no external torques, parallel to the z 
axis, acting on the system. 

These equations are not difficult to model numerically and can be solved by standard techniques. 
Although a brief analysis of the above equations shows that if the particles are travelling at speed 
greater than the escape velocity of the system, then the initially vertical path of the projectile will 
quickly turn into a "horizontal" path across the face of the accretion disk. A computer simulation of 
this phenomenon is given in Fig. 9 ( see also Liffman and Brown (1995)). 

In Fig. 9, particles are ejected at r = 0.04 AU with different vertical velocities ranging from 149 to 
174 km s^ 1 . The particles are assumed to move out of the gas outflow at r = 0.1 AU, whereupon they 
encounter the gas halo of the accretion disk and their subsequent motion is governed by gas drag plus 
the gravitational force from the protostar. 

As the particles move through the halo gas of the accretion disk, they will acquire, by gas drag, 
angular momentum from the disk. If we assume that the centrifugal force of the halo gas balances the 
radial component of the protostar's gravity, then the angular speed of the halo gas, vg gas , will be given 
by 

_ rVGM 

Ve,gas ~ ^ 2 + ^2)3/4" ^ i6 > 

Once the particle has come to rest, relative to the halo gas, it will have the angular velocity given by 
Eq.(73) and a specific angular momentum, h, given by h = rvg jgas . As can be deduced from Eq.(69), 
such a specific angular momentum implies that r = 0, and the only force acting on the particle will be 
the z component of the gravitational force, which will point towards the accretion disk. As the particle 
moves towards the accretion disk, the gas density and angular velocity of the gas will increase, thereby 
keeping r w 0. As can be seen from Fig. 9, the subsequent path of the particle is roughly parallel to 
the z axis. 

The paths of the captured particles, shown in Fig. 9, can be approximated to that shown in Fig. 
10, where the ascending path length, I , is given by I = Kx, with K being a number in the range of 
4 ± 1, (Liffman and Brown 1995) and X is the "stopping distance" as defined by the equation, 

X = ^EL. (74) 

As a justification for the above equation, we note that a macroscopic particle will come to rest when 
it has encountered a total gas mass approximately equal to its own mass. The "stopping distance" , so 
defined by this prescription, is easily shown to be the X length scale as given in the above equation. So, 
for a constant gas density, the larger or more dense a particle is, the further it will be able to travel 
before it comes to rest. 

The range of a recaptured projectile is simply given by the formula 

r^lcos(9), (75) 

where 9 is the angle between the ascending path of the particle and the midplane of the accretion disk. 

The simplicity of Eq.(75) suggests that aerodynamic size sorting may occur to particles that are 
recaptured by the protostellar system. To see this, we suppose that we have two particles with different 
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Particle Trajectories with Gas Drag 
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Fig. 9. — Simulation data showing chondrule-like particles travelling above the disk and falling back 
into it. A solar-mass protostar is located at r = 0, z = 0. Surrounding the protostar is an accretion 
disk, the scale height of which is shown in profile. Chondrule-like particles (radius = 0.1 cm, density = 
3.5 g cm -3 ) are given a velocity boost in the z direction, the magnitudes of which (in km s -1 ) are shown 
next to the trajectories. Particles that are subject to sufficiently high gas drag are later recaptured. 




Fig. 10. — An approximation of chondrule trajectories as a triangular path. 
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mass densities: p\ and p%. Suppose further that these two particles fall back to the accretion disk at 
the same distance from the protostar. Then we can write the radius ratio of particle 1 to particle 2 as 

^ - ^Q(r) , (76) 

«2 Pi 

where 

Q(r) = ^i^i^M (77) 
Q[r) K, p g2 cos(^) ■ {<<) 

The factor Q(r) is dependent on the initial ejection speeds of the particles (through K and 9), the 
initial distance of the particles from the protostar ( again K and 9) and the scale height of the accretion 
disk (through p g ). If the particles are created at about the same distance from the protostar then we 
would expect that K, 9 and p g would be similar for both particles, since the average flight path of both 
particles would be about the same. This would imply that Q(r) « I. So, from the collection of particles 
that fall to the accretion disk at r, the denser particles should have smaller radii. 

As has been discussed, size sorting is density dependent in that dense Fe-Ni chondrules are always 
smaller than less dense silicate chondrules. The mass densities of the two types of chondrules are (Skinner 
and Leenhouts 1993) psi ~ 3.8 g cm -3 , and pp e ~ 7.8 g cm -3 , which implies that p,si/pFe ~ 0.5, and 
so if Eq.(69) is applicable to chondrules, we should expect that 

— « 0.5Q(r) . (78) 

asi 

Using published data for the meteorite types H, L, LL, (Dodd, 1976) and CR (Skinner and Leen- 
houts 1993), we can compute the average Fe to Si size ratio for all these different types of meteorites 
(excluding Bjurbole) 

<QFe> \ =0.52 ±0.16. (79) 

< a Si > I H,L,LL,CR 

The corresponding approximate mean Q values is 

<Qh,l,ll,cr>= 1-04 ±0.32 . (80) 

Of course, the agreement between the above Q value and our theoretical model should be treated 
with caution, since we have only presented the bare beginnings of a quantitative model. Nonetheless, 
it does illustrate the potential of the "Jet" model to explain the phenomenon of size-sorting. 

Finally, we return to our discussion of angular momentum transfer and disk accretion. Let us 
consider a ring of material in an accretion disk at a distance r from a protostar. The angular momentum 
of this material is 

L(r) = m(r)VGMr , (81) 

where m(r) is the mass of the ring of material at r. Now suppose that material falls onto the accretion 
disk, and this infalling material has essentially zero angular momentum. In such a case, the angular 
momentum of the ring is conserved, and the mass of the ring becomes a function of time, i.e., m = m(r, t). 
For such a case, we can differentiate Eq. (81) to obtain 

dr f^ r \ d m 



> , • (82) 
dt \mJ dt v 1 
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or 



r 




(83) 



Thus, if infalling material doubles the mass of the ring, it will move from its initial position r(0) to 
r(0)/4. It is via this mechanism of mass and angular momentum transfer that disk accretion may be, 
in part, mediated. Indeed, because it is the halo gas of the accretion disk that will be transferring most 
of the angular momentum to the "Jet projectiles", it is possible that the upper layers of the accretion 
disk are the ones that undergo most of the accretion, leaving the midplane relatively untouched. 

This type of process may explain an implicit contradiction between observations and meteoritics. 
Observations suggest that protostars keep on accreting material from their disk for periods of up to 10 7 
years (Cabrit et al. 1990). Radiometric data from the decay of 129 I and 26 Mg suggest that meteorites 
accreted material for periods of order 10 6 - 10 7 years. How could the meteoritic material have been 
preserved if a major portion of the solar nebula was accreted onto the protosun? The answer, we suggest, 
is that disk accretion was, in part, altitude dependent. Material flung from the protosolar jet mediated 
the angular momentum transfer and one component of this mass transfer was the chondrule. 



In this paper, we claim that chondrules have the physical characteristics expected of ablation 
droplets that have been formed and ejected by a protosolar jet, and then recaptured by the solar nebula 
through the action of gas drag. 

Such a model, however, is immediately confronted with the complementary composition of matrix 
and non-matrix material in meteorites. For example, Wood (1985) discusses the case of Murchison, 
where the matrix has an Fe/Si ratio of 1.23, while the non-matrix material (chondrules, CAIs, isolated 
crystals) is ~ 0.2. These two dissimilar components combine to give an Fe/Si ratio of 0.81, which is close 
to the solar value of 0.9. This is unlikely to be accidental, and is clear evidence for the local formation 
of chondrules, and the refractory component of chondritic material. 

Before we throw out our wind-transport model, however, one should note that chondrules, at least 
in this model, form in a relatively small region of the solar nebula, i.e., in or near the boundary layer 
of the protostellar system. If the chondrules have a low Fe/Si value then the surrounding material will, 
by mass balance, have a high Fe/Si ratio. If this latter material comes in a non-gaseous form ( e.g 10 
fim dust, CAI &c) then it too will be ejected with the chondrules. One will obtain the desired mass 
balance if all this material lands back into the solar nebula, in a uniform manner, and over a long period 
of time. 

For such a model to work, we require a fairly large component of matrix material to be made from 
dust that has been recycled through the protostellar jet. There are at least two consequences if this 
idea is correct. First, this recycled dust would have to be more refractory than CI material, since the 
jet formation region of the inner solar nebula would have been far warmer than the regions where most 
chondrites were formed. Second, if matrix material were formed from dust that had been lofted into 
the upper atmosphere of the solar nebula then this may be an observable phenomenon. 

To understand this latter point, we need an estimate for the amount of dust that should be resident 
in the upper atmosphere of the solar nebula at any particular time. 



6. 



CHONDRULE-MATRIX COMPLEMENTARITY 
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Protostellar jets have average mass loss rates of order 10~ 8 M yr _1 and since the mass of dust 
to gas in the Interstellar Medium (ISM) is 1/100, this implies that <, 10~ 10 Mqjt" 1 of dust is blown 
out by the protostellar jet. So, an upper limit for the amount of dust that is lobbed into the upper 
atmosphere of an accretion disk in a Classical T Tauri Star (CTTS) per year is 10~ 10 M Q . 

The dust settling timescale (r sett j e ) can be deduced from Eq. (15) with v gz = 0, i.e. , 

C D . 2 GMrripZ 

m P Z = ~^P9 Z A P ~ p + ^3/2 • ( 84 ) 

The motion of the dust is subsonic, so the drag coefficient takes the Epstein form for Co (Eq. 18). To 
compute the velocity of the particle, we note that z«0, which implies 

Pg v 

where uk(t,z) = V GM/[r 2 + z 2 } 3 ^ 4 is the Keplerian angular velocity at the point (r,z), and v = 
^8kTg/nm g is the mean Maxwellian speed of the gas particles. 

If z <C r then uk(t, z) ~ wjr(r) = a/ GM/r 3 and we can compute the settling time: 

r, Mle « 8, 000 (ft/10- 11 gcm-'y/l l-m'-'^/l ™f yr (M) 

(ap/l^m)(pp/l g cm d ) 

Thus, we have characteristic dust-settling timescales in the range 10 3 — 10 4 years, which means that up 
to 10~ 7 — 10~ 6 M of dust will be in the upper atmosphere of a CTTS accretion disk at any one time. 

As discussed in Natta (1993), 10~ 7 M of high-altitude dust may produce the observed "flat- 
temperature distributions" in CTTSs. These temperature distributions are a surprisingly common 
phenomenon in CTTSs. They arise when the temperature of the disk does not decrease as rapidly with 
distance from the protostar as one would predict from standard accretion disk theory. Natta suggested 
that a spherical halo of dust around a CTTS would reflect light from the protostar and into the accretion 
disk, thereby increasing the temperature of the outer disk. 

Of course, the idea that protostellar jets can loft dust into the outer parts of the surrounding 
accretion disk requires a quantitative investigation to determine whether dust grains can be ejected 
from the jet flow. There will be a size limit where particles smaller than a certain size will simply be 
entrained in the jet flow and ejected from the system. We simply note, that we require dust ejection 
similar to that shown in §5 to account for chondrule-matrix complementarity. 



7. CONCLUSIONS 

Chondritic meteorites are typically an agglomeration of igneous rocks, i.e., chondrules and refractory 
inclusions (formation temperatures 1500-2000K) , surrounded by sedimentary material that, in some 
cases, has never experienced temperatures greater than 500K. This unusual structure has prompted 
theorists to develop heating mechanisms (e.g. lightning) that can provide brief, intense impulses of 
energy in the otherwise cold outer regions of the solar nebula. These energy impulses are presumed 
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to have melted small dust aggregates into chondrules, which were then incorporated into larger dust 
aggregates that eventually formed meteorites. 

Such energy impulse theories are not required in the 'Jet' model of chondrule formation. Chondrules 
are formed in the hot inner regions of the accretion disk adjacent to the protostar. There is no difficulty 
in obtaining the required temperatures, because in or near this region the accretion disk dumps around 
half of its gravitational energy. The protostellar jet ejects the chondrules from the hot inner disk and 
gas drag brings these particles back to the cooler outer regions, where they could be incorporated into 
growing aggregates of cool nebular material. The "plum pudding" structure of chondritic meteorites is 
a natural consequence of this model. 

One of the many problems in chondrule formation is the deduced low cooling rate (~ 1 - 1000 
K/hour) for these particles. This cooling rate is many orders of magnitude smaller than that expected 
for a particle radiating directly into space. In our model, chondrules are formed in the optically thick 
regions of the inner accretion disk. They will, therefore, be in thermodynamic equilibrium with the gas 
and their temperature variations will be damped. 

Of course if the chondrules were to remain in this environment, their temperatures would not 
decrease. Chondrules, however, are produced by the ablative interaction between a streaming gas flow 
(perhaps produced from the merging of magnetic field lines) and molten material. Particles, that are 
small enough, will be swept up with the flow. As the gas flow moves away from the midplane of the 
inner accretion disk, it will expand and therefore probably cool. Particles that are entrained in this gas 
flow will also cool at the same rate. 

Particles that are moving with the gas flow may be ejected from the accretion disk, but it is 
possible that a particle may simply hover at some distance away from the midplane of the disk. If 
the gas density in the flow is below a critical gas density (~ 10~ n g cm -3 ) then the particles will 
undergo damped oscillations around their hover points. It is due to these oscillations, we suggest, that 
chondrules can undergo reheating. Smaller particles of the same mass density will still move past the 
hovering particles, and so may collide with these particles. 

Such a scenario allows for the formation of adhering compound chondrules. These chondrules pose 
a major problem for chondrule formation theories as it is nearly always the smaller chondrules that 
were plastic at the time of collision. Smaller particles will lose heat more readily than larger particles, 
since they have a larger surface to area ratio. If chondrules were formed in a single flash heating event, 
we should expect the smaller particles to become solid before the larger particles. That we actually see 
the opposite behaviour (i.e., the large particles were solid, while the small particles were semi-molten), 
strongly suggests that single flash-heating models require some modification. 

In the Jet model, adhering compound chondrules are formed when small chondrules collide with 
larger chondrules that are stationary in the flow. The smaller particles that are moving with the flow will 
be warmer than the larger hovering particles, because the latter particles have had time to equilibrate 
with the local temperature of the gas. If the temperature gradient in the flow is sufficiently steep (~ 10 
to 100 K/km) at this collision point, compound chondrules with the observed structure will be formed. 

Once a compound chondrule has been formed, and has stopped spinning, the smaller secondary 
chondrule will orient the entire compound chondrule, just like a weathervane, such that the secondary 
chondrule is pointing in the direction of the flow. Any incoming secondary chondrule will hit the rear 
end of the primary chondrule. The resulting triple compound chondrule, will have two secondaries 
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chondrules that will tend to "avoid" each other. Of course, this prediction implicitly assumes a faily 
steady gas flow, an assumption that may be incorrect. However, if this prediction is found to be valid, 
it would be major piece of evidence in favour of the Jet model, as it is difficult to produce such an effect 
with other chondrule formation theories. 

The fragmentation of chondrules is a simple extension of compound chondrule formation, where 
instead of one particle being plastic at the time of collision, we now have two solid particles colliding at 
higher velocities. We expect the fragmentation zone to be 'above' the compound-chondrule formation 
zone, since the gas flow will be cooler and the flow speed higher as one moves further away from the 
midplane of the disk. 

The chondrule size range is yet another aspect of chondrule formation which has not been satis- 
factorily explained. In the Jet model, the sizes of the ejected droplets are determined by the balance 
between gas drag and the surface tension. A jet flow with high energy density will tend to make molten 
droplets smaller, while molten materials with higher surface tensions will tend to form larger droplets. 
Observations suggest that protostellar jets are produced within 0.1 AU of solar- mass protostars. The 
minimum energy density of a wind that can eject Si and Fe droplets from such close proximity to the 
protostar is such that the radii of these droplets is < 1 cm. Protostellar winds that form larger droplets, 
will have a lower energy density and will not be able to eject them from the jet formation region and 
we will not see them in meteorites. 

This ejection mechanism may also explain the complementary chemical structure of chondrules and 
their surrounding matrix. For this to occur, we require that refractory dust as well as chondrules are 
ejected by the jet flow so that mass balance will be obtained once the dust settles back to the solar 
nebula. This leads to a simple calculation, which suggests that up to 10~ 7 — 10~ 6 M of dust will be in 
the upper atmosphere of a CTTS accretion disk at any one time. These figures appear to be consistent 
with those deduced from observation (Natta 1993). 

The transfer of dust and chondrules from the inner to outer portions of the accretion disk, necessarily 
requires the removal of angular momentum from the upper atmosphere of the accretion disk. This 
will increase the viscosity of the disk and hasten the stratified infall of disk material, i.e., the upper 
atmosphere of the disk will accrete more readily than the central plane of the disk. 

Finally, aerodynamic size sorting of the ejected particles is a natural consequence of the Jet model. 
Ejected particles that are captured by the accretion disk must suffer aerodynamic size sorting, since the 
particles will be subject to gas drag. Particles that are not subject to gas drag will simply leave the 
protostellar system, since their initial velocities were higher than the escape velocity of the system. 

It is for all these reasons that we consider chondrules to be ablation droplets formed by a protosolar 
jet in the first 10 6 — 10 7 years of the solar system. 
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Jet Flow 




Fig. 11. — A pictorial summary of the ideas presented in this paper. An accretion disk surrounds a 
protostar. At a distance within 0.1 AU of the protostar, the disk produces a protostellar jet (only one 
side of which is shown). Infalling km-sized bodies move into this hot "jet zone" and are subsequently 
ablated by the jet wind. The resulting droplets, if small enough, move with the flow and away from 
the ablating body. It is our claim that these ablation melt-droplets eventually become chondrules. 
At higher altitudes, within the accretion disk, some of the larger particles may hover in the gas flow. 
Smaller particles will still move with the flow, however, and may collide with the hover particles. 
These chondrule-chondrule collisions produce compound chondrules and chondrule fragments. The 
fragmentation of chondrules occurs at higher altitudes relative to the compound- chondrule formation 
zone, because the flow cools as it increases in altitude and the two colliding chondrules will be solid 
at the time of collision. If the protostellar jet is sufficiently powerful, it will eject particles from the 
inner accretion disk, and if the orbital angular momentum of the ejected particles is large enough, the 
particles will move out of the jet flow and travel across the face of the accretion disk. With sufficient gas 
drag from the upper atmosphere of the accretion disk, these (now size-sorted) particles will be brought 
down to the outer accretion disk, where they will be incorporated into planetesimals. 
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